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Abstract 

By the method of invariant manifold, we investigate the Ito equa- 
tion numerically with high precision. By the numerical results, we can 
completely determine the form of analytic soliton solutions for the Ito 
equation. In fact, by the numerical data we have succeeded in deciding 
the analytic form of the r-function, which is more general than the old 
assumptions. This may suggest that we should think more deeply about 
what the soliton is. 
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1 Introduction 

The Ito equation [1] was first introduced in 1980 as the Hirota bilinear 
equation 

D t (D t + D 3 x )f-f = 0. (1) 

By introducing variable u = 2(\nf) xx , Equation ([T]) is transformed to a soliton 
equation 

u t t + u xxxt + 3(2u x u t + uu xt ) + 3u xx / u t dx = 0. (2) 



Non-local equation ^ is usually written into a local one by the transforma- 
tion u(x,t) = u(x, —t), -^u(x,t) — 3-^v(x,t). Neglecting the tilde symbol and 
the integral constant, we get 

u t = 3v x , 

v t = (uv) x + v xxx . (3) 
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In the sequel system ([3J will be referred to as Ito equation. In Ref.[2], Drinfeld 
and Sokolov also found the equation (|3j in their study on the relationship be- 
tween soliton equation of KdV type and Lie algebras, thus some authors referred 
([3J as Drinfeld- Sokolov system [3J. The Hamiltonian structures of Ito equation 
was obtained by Liu [3]. 

The Lax pair for © is [3] 

L = d 3 + ud + u x + vd- 1 , 

P = d{d 2 + u). (4) 
Also in PQ Ito himself gave a solution of (p} 

/(N) N \ 

f = H exp I H ^tf MiMj + (MVi i (5) 

M=0,1 \i<3 i=l / 

where 

?7i = + flji + 77°, 
«i = -P?, 

fa -Pj)(p? ~P 3 3 ) 
(Pi+Pj)(j>i +P 



expCy = 7— — — T7 -3— 3-. (6) 



It had been a long time taking ([5]) as the r-function of the soliton solution 
of ©. But one has good reason to question whether ([5]) covers all soliton r- 
function of the Ito equation: literature \Y\ has assumed that for soliton solutions 
/ is a sum of exponential traveling waves. Zhang and Chen [5] are among the 
first ones to try to enlarge the class of soliton solutions of the Ito equation. It 
seems to us that the solutions obtained in [S] are still too limited. Li et al. [6 
had investigated the Ito equation by the method of invariant manifold (IM) and 
solved the first two constraints completely. But they [B] did not give the general 
solution for the higher order constraints. The results of this paper will reveal 
some important structures for higher order constraints discussed in [6]. 

This paper is organized as follows. Section 2 will derive the ordinary dif- 
ferential equations (ODEs) for solving the IMs suggested by [6]. The ODEs in 
Ref.[B] is not suitable for numerical integration. Here we give the suitable one. 
Also here we add an equation for f , which is crucial to analyze the IMs. Section 
3 will present the numerical results for the ODEs obtained in section 2. These 
numerical results will reveal the general structure of r function corresponding 
to the 2-soliton and 3-soliton solution of the Ito equation. Section 4 is the 
conclusions. 

2 ODEs for the time evolution of solitons of the 
Ito equation 

For brevity we use the following notations: /' = fj) / = §{• The definition of 
d- 1 is d^d = dd- 1 = 1 0. 
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With eigenfunction ipi, the Lax pair becomes Ltpi — Aj^i an d tpi — Pipi, 
which are equivalent to the following differential equations 

i>i = + (u^y. (7) 

According to [6], the constraint for Equation (J7|) is 

n 

v = ^aiipi, (8) 
i=l 

where djS are real or complex numbers and we can take ai = 1 without loss of 
generality. 

Let ipi = f'i- Then (JTJ) becomes 



Ai^i - p'I" - ( u p'i)' = v (Pi + h{t)), 

pi = p'-' + Wi + Ci(t). (9) 
= and Cj(t) = 0. 



It is easy to prove ki(t) 
Altogether we have 



hp\ ~ P?' ~ = vpi, 

Pi = Pi' + UPi, 

71 

i=l 

u = 3v'. 



(10) 

(11) 

(12) 
(13) 



Remark 1: From Equations (fit))) and flT2j) , we can easily see that with re- 
spect to x there are n equations with n+1 unknown functions: {ipi, ip2, • ■ ■ , <£>n, u} 
So we can select u(x, i)|t=o = uq(x) to be an arbitrary function, i.e., with respect 
to x, u is not limited to the form 2(ln^ a^e iX ) xx DP- 

By ([TP]). (|TT|) . dHJ) and ([15]). the time evolution of u, tfi, ip", <p"' become 
a closed system of ODEs: 

u = 2>v' . 

Pi = Pi' + u P'i- 
& = (tf' + uip'J 

= \iip' t - Vlfii, 

Pi = ( X iPi - v Pi)' 

= hp" ~ v'pi - Vip'i, 

<£' = {Xip'l - v'tfi - vtpi)' 

= A^f - v"cpi - 2v'i P ' i - vcp'f, (14) 
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where i = 1, 2, 3, • • • , n. Clearly the ODE system (fT4")l contains 4rt + 1 equations. 
Next we will introduce an extra differential equation 

n 

f = X>. (15) 

4 = 1 

Equation (| 15|) is very important to reveal the solution structure of Equations 
Summing (fTOj) over i, we immediately get 



Integrating (1161) we get 

/ „ \ 2 



5lX> +E(A^-^"-^)- (17) 



2 , 

\i=i / »=i 



It is easy to prove 7 is independent on t. We will use the conservation of 7 to 
illustrate the precision of the numerical solutions. 



3 Numerical results 

Because the case n = 1 had been completely solved in [B], we start numerical 
study from n = 2. 



3.1 n=2 

The parameters A 2 , i = 1,2, can be taken any value in Equations (|14l) . For 
simplicity we take 

A x = l, A 2 =2. 

The initial values for unknown functions in (fl4l and (|T5|) are also arbitrary. 
In order to make the numerical solving procedure simple we randomly choose 
them as 



Vi (0) = 


1. 


¥> 2 (0) = 1, 


ri(o) = 


1, 


^(0) = -1 


^i'(O) = 


1, 


¥>a(0) = 1, 


<(0) = 




^"(0) = 2, 


u(0) = - 


-3, 


f (0) - 0. 



(18) 
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Fig. 1 The difference between j(t) and 7(0). 

The calculated 7 is very close to 7(0) = —3. Fig. 1 shows that the numerical 
precision in our numerical integration for the ODEs (|14p is very high. 

Numerical results show that for lots of initials tpiS are asymptotically to 
constants. For initial values (fT8l). the numerical results are 



ifl(+oo) 
tp 2 (+co) 



(100) 

^ 2 (100) 



-8.515821092114916270761375823327277321106, 
15.71501739387443269240781890998386797269. 

(19) 
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of ip[ and tp' 2 are determined by 



M=( A ^-^(+f -^+~) ). (20) 



The eigenvalues of — M are 



O x = 3.1444367053092452142565168204216716495466, 

fi 2 = 1.0547595964502712073899262662349190020374. (21) 



Define 



uix = -Qi, Lu 2 = -0,2. (22) 
By the numerical values of (pi(t) and f2(t), t £ (0, 3), we can determine 

/,s _ i,je{i -1} ' , s _ -1} ' , , 

ije{i,-i} ' i,je{i,-i} 



where we can assume 6_i,_i = 1. Note that ()23[) are obtained by only a fuzzy 
thinking. The calculated a}j, a 2 • and 6ij are 

a ia = -0.66605554882019023429631637563968997173, 
4 _! = 0.27052090248727273290847449418182448488, 
ai^ = 3.6002398946839792257416422206042247363, 
a 1 _ 1 _ 1 = -0.117428488595883427468489650014095681785, 
a? ;1 = 1.229132742665024456496131959076675712559, 
a\ _ x = 5.796357753327621304524412601597533898065, 
a\ x = -2.856445923073834469957280434899450711875, 
al x _! = -1.081767813163632994177953436642495030883, 
6i,i = 0.078213896418856586401163070522148337413, 
6i _i = 1.19199676640685627533532926755144489269, 
= 0.81706609692946543514881835105867033706, 
6-1,-1 = 1. (24) 
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Fig. 3 The difference between tpi and (p\. (pi denote ipi with ([241 
and (p\ denote the numerical solution of tp\. 



Fig. 3 shows that the form of tp\ must be (|23|) . 

By the numerical values of f (t), t 6 (0, 3), we can determine 



f(t) = c + cii + In ( £ 6 i , j e iu,lt+iW2t 
,;,je{i,-i} 



and the values of cq and c\ are 



c = -2.254578789614522148435364381123302299, 

ci = 3.000000000000000000000000000000000005 w Ai + A 2 . 
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Fig. 4 The difference between / and /. / denote f in (f2"5"j) with (f2Ti|) 
and / denote the numerical solution of f . 



Fig. 4 shows that the form of f must be ([25 
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3.2 n=3 



The procedures for n = 3 is just the same as the one for n — 2. All 
numerical errors are also far smaller than 10~ 30 . We will omit the plots here. 
We arbitrarily take the parameters Aj, i — 1, 2, 3 as 

Ai = 1, A2 = 2, A3 = 5. 

Then we set the following initial values 

¥>i(0) = l, <p 2 (0) = 1, ¥> 3 (0) = 1, 
^(0) = 1, ^ 2 (0) = -1, ^(0) = -3, 

(0) = -1, ^'(0) = -l, ^'(0) = -l, 
<(0) = -2, ^"(0) = -2, ^"(0) = -2, 
u(0) = -6, /(0) = 0. (27) 

We can obtain 

931(00) « ^i(72) = 14.92110270021879759145408479891323291706, 
932(00) « ^ 2 (72) = -41.37915928188619028119192656047612164361, 
933(00) « 933(72) =44.53116386683746280518744342554027669437. (28) 



Ai - 931(00) -931(00) -931(00) 
M= I -932(00) A 2 - 932(00) -932(00) ). (29) 
-933(00) -9^3(00) A 3 -933(00) 

The eigenvalues of — M are 

fii = 6.277038684404317215906811766255007511305, 
ft 2 = 2.341721658570986197790603783195012345933, 
ft 3 = 1.454346942194766701752186114527368110582. (30) 

Define 

wi = ^fii, lo 2 = ^f2 2 , w 3 = 7^3- (31) 
By the numerical values of <pi(t),t E (0,3), we can determine 

i,j,fce{i,-i} 



8 



where we can assume = 1. The calculated a\ - k and hj^ are 

a\ xl = 0.529262661391713354983735535897686104404, 
a\ tl _! = -0.350348759673940084264112387560739906945, 
a\_ 1A = -8.2675670367386623516544803221517811099, 

= -0.0254209498663860215070864228788457306588, 
a 1 :!,!,! = -8.2069634233611868635668643982404949776, 
_i = 0.94894068738511728728577621078955537731, 
a?_i _ M = 25.8865938214425998112477469386616735198, 
ai-i _i _i = -0.80423010558221252574876705199352366617, 
6i,i,i = 0.035470747170981702373487982709002505528, 
6i,i -i = 0.126837598026745943988081084806178876739, 
6i = 1.38001767415062124882858267790796406951, 
61-1,-1 = -0.0229217114935723223611012392981108059801, 
fo_i,i,i = 0.758483231785985808547518495835409278743, 
6-1,1,-1 = 0.473751827832812619094405400960157166769, 
6-1,-1,1 = 5.95862752752346760630497369960292851901, 
6-1,-1,-1 = 1. (33) 

Similarly the other of • k and a? ■ fc ( the coefficients of <^2(i) and <£3(t) ) can 
also be calculated. But they do not have much importance, and we will not list 
them here. 

By the numerical values of f(t), t G (0, 3), we can determine 

f (i) = c + at + In I J2 b id , k e iuit+ju > t+ku * t \ , (34) 
\i,j,fee{i,-i} / 

and the values of c and a are 

c = -4.5463675370738389259379741125305027486, 

a =8 = A!+A 2 + A 3 . (35) 

4 Conclusions 

By the numerical results in Section 3, we can see that for general n 

2 
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It can be shown that cq and 6» x » 2) ... ,i„ are all dependent on x and that u>i 
i = 1, • • • , n are not dependent on x and t. Therefore we can write 




2 

h^ 2 ,..^{x)e^ t+ ^- t+ - +M ) , (37) 



where the demand 6_x,— l,--- l = 1 has been removed. By (|15p we can only get 
fzt = v. The r-function of the Ito equation satisfies r xt = v and t xx = ^u. So 
we can not immediately identify r with f . But we can still show r has the same 
form as (1371) or we can write 



r = ln & ll , l2 ,-, 4 „(a;)e ll " lt+42tJ2 * + - +t """ t . (38) 

\ii,ta,— >*n6{l,-l} / 

Of course, 6j li j 2i ... t i n (x), i = 1, ••• ,n, in (|38p are not arbitrary. They must 
satisfy some constraints. In fact we know there is only one arbitrary function 
among all the ^ 1 ,i 2 ,..-,i„(2;)s. 

From the numerical results in Section 3, we can also see that the eigenfunc- 
tion ipj must take the form 

y~* q3 '• ^ a ,^gii"it+i2"2*H H,i"„t 

Vj = t " t2 '-''" g{1 ^ 1} -' 1 ' t2 ' . (39) 

J ^2 bi i i (a;)e iltJl * _ '"* 2 ' i ' 2t ^ — \-ini^n,t v ' 

iiifer" ,inG{l, — 1} 

A discovering of this paper is that the ODE system (| 14|) is a completely 
integrable one, which seems to be never seriously studied in the literatures. 
Therefore we can conclude that the time evolution for the multi-solitons of the 
Ito equation is integrable. The dependence on x for the multi-solitons of the 
Ito equation is Equations (|10|) and (fT2|) . which are not completely solved yet. 
Recall that the KdV equation with vanishing boundary condition is completely 
solved by the multi-soliton solutions of the KdV equation. Likely we conjecture 
that the Ito equation with vanishing boundary condition is completely solved 
by the multi-soliton solutions governed by Equations (flOl) . (|TT|) . (fl2|) and (|13j) . 

The method proposed here provide an interesting way from numerical studies 
to analytical ones. We expect this method can also be applied to some other 
nonlinear systems. 
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